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Abstract 

Chamfer distances play an important role in the theory of distance transforms. 
Though the determination of the exact Euclidean distance transform is also a well 
investigated area, the classical chamfering method based upon "small" neighbor- 
hoods still outperforms it e.g. in terms of computation time. In this paper we deter- 
mine the best possible maximum relative error of chamfer distances under various 
boundary conditions. In each case some best approximating sequences are explicitly 
given. Further, because of possible practical interest, we give all best approximating 
sequences in case of small (i.e. 5x5 and 7x7) neighborhoods. 
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1 Introduction 



Suppose we measure distances between grid points of a two-dimensional grid 
and we want to approximate the Euclidean distance by a distance function 
which can be computed quickly, without calculating square roots. We may then 
use the class of chamfer distances. They are obtained by prescribing the lengths 
of the grid vectors in a so-called mask M p := {(x,y) G Z 2 : max(|a;|, \y\) < p} 
(for some positive integer p) such that the values at (±x, ±y) and (±y, ±x) 
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are all the same, and by defining the length function W as follows: the length 
W(v) of any vector v G Z 2 is defined as the minimal sum of the lengths of those 
vectors from M p , repetitions permitted, which have sum v. The literature on 
chamfer distances is very rich. See Borgefors [Tf2f3] for the basics, [4"f5] for lists 
of (2p + 1) x (2p + 1) neighborhoods for 1 < p < 10, and [6] for an overview 
of applications. Further, recently many related results have been obtained by 
several authors, concerning distance transforms and their explicit calculation 
using different kinds of neighborhoods in certain (mostly 3D) grids. For ex- 
ample, Strand, Nagy, Fouard and Borgefors [7] gave a sequential algorithm for 
computing the distance map using distances based on neighborhood sequences 
in the 2D square grid, and 3D cubic and so-called FCC and BCC cubic grids, 
respectively. Similar results for other kinds of grids are also known, see e.g. [8] 
(rzD hexagonal grids), [9] (diamond grid) and [10] (general point grids) and 
the references given there. 

Classical chamfer distances using 3 x 3, 5 x 5 and 7x7 neighborhoods given 
by Borgefors [Tf2] are generated by the masks 

14 11 10 11 14 
4 3 4 11 7 5 7 11 

3 3 , 10 5 5 10 and 

4 3 4 11 7 5 7 11 

14 11 10 11 14 



51 43 38 36 38 43 51 
43 34 27 24 27 34 43 
38 27 17 12 17 27 38 
36 24 12 12 24 36 
38 27 17 12 17 27 38 
43 34 27 24 27 34 43 
51 43 38 36 38 43 51 



respectively (with the actual generator entries underlined). For comparison 
with the Euclidean distance the values of the neighborhoods have to be di- 
vided by 3, 5 and 12, respectively. The approximations to a/2 m 1.41 are 
therefore 4/3 ~ 1.33, 7/5 = 1.4 and 17/12 « 1.41, respectively. For alter- 
native neighborhood values see Verwer [TTlfT2] . Thiel [13], Coquin and Bolon 
[H], Butt and Maragos [15] and Scholtus [4j. More specifically, in [2] the 
minimization of the error between the Euclidean distance and the local dis- 
tance was considered over circular trajectories similarly to [TT|T2] rather than 
linear ones [2lTl~3~] . The approximation error can also be measured based on 
it is done in [15] with calculating the difference between a disk of large 
size obtained by chamfer metric and a Euclidean disk of the same radius. The 
determination of the exact Euclidean distance transform is also a well investi- 
gated area (see e.g. [T6ll6lll7fl8|rT9] ). but the classical 3x3 chamfering method 
still outperforms it in terms of computation time and simple extendability to 
other grids. 
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In this paper we determine chamfer distances best approximating the Eu- 
clidean distance in a certain sense. In each neighborhood size some best ap- 
proximating sequences are explicitly given. Further, because of possible prac- 
tical interest, we give all best approximating sequences in case of small (i.e. 
5x5 and 7x7) neighborhoods. 

Throughout the paper, as a measure for the quality of a length function W 
defined on Z 2 we use the so-called maximum relative error (m.r. error for short) 



E := limsup 

\v\—>oo 

where |.| denotes the Euclidean length. The Mi-, M2- and M3-neighborhoods 
given above yield rounded _E-values 0.0572, 0.0198 and 0.0138, respectively. 
Firstly we shall prove that the smallest possible constant E p for the mask M p 
under the condition that W(x, 0) = \x\ for x G Z is given by 



% = 2 = 2 rUl-r]. 

' P Z p Z \P J 

In particular, Ef « 0.0551, E$ « 0.0187 and Ef 0.0089. Comparing these 
values with the E- values given above, one can see that the E^- values yield 
approximately 4%, 6% and 35% improvement, respectively. The B refers to 
Borgefors who was the first to consider such neighborhoods. 

Secondly we consider the case D in which W(v) > \v\ for all v G Z 2 . (The D 
refers to the fact that W(v) dominates \v\.) The optimal m.r. error under this 
restriction equals 



In particular, E? w 0.0824, E% w 0.0275 and Eg w 0.0131. 

Thirdly we shall prove that the optimal E'-value without any restriction on 
the neighborhood defined on M p (i.e. dropping the condition W(x,0) = \x\ 
for x G Z) equals 



c 



2p 2 + 2- 2py/p r TT - 1 



2p 2 + 2 - 2pvV+T + 1 




In particular, w 0.0396, E^ 7 w 0.0136 and £f 0.0065. In 1991, on using 
the symmetry in case C the value of E^ was computed by Verwer [TTyT2] in 
terms of trigonometric functions. The C refers to the word central. In 1998, 
because of geometric considerations, Butt and Maragos [15] chose to use the 
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error function 



W(v) 

which of course is small if and only if is small. In general it gives different 
error values, but the values for E& are equal to the values obtained by the 
above error function (cf. Scholtus [1]). We prove the correctness of the above 
Ep values. In doing so, our motivation is twofold: on the one hand, by a sim- 
ple reasoning we obtain these values immediately from the values of E^, and 
on the other hand, our proofs are mathematically rigorous while the corre- 
sponding arguments of Verwer and Butt and Maragos contain some hidden 
assumptions. Namely, by certain plausible but not explicitly verified geometric 
arguments they restrict their attention and investigations to certain values of 
the neighborhoods in question, and they perform exact investigations only for 
these values. 

We shall further study an auxiliary class of neighborhoods on M p , viz. the 
class of neighborhoods satisfying N c (v) = oo for all v = (x,y) G M p with 
either x < p or y < 0, M c {v) — p for v — (p, 0), and M c {v) = c\v\ for v = (p, k) 
with < k < p. Here c is a constant close to and at most equal to 1. Informally 
speaking, the use of such neighborhoods means that only such steps (t>i,t>2) 
are allowed, where v\ is a positive multiple of p and v-i is nonnegative. Further, 
beside Af c (p, 0) = p the weights of the other such neighborhood vectors are 
their Euclidean lengths, multiplied by a factor c < 1. All the other vectors 
of the neighborhood are forbidden to use, thus they have weights oo. For 
example, the weights for the neighborhood J\f c with p = 2 (i.e. for M2) are 
given by 

00 00 00 00 c\/8 
00 00 00 00 cv^5 
00 00 00 00 2 
00 00 00 00 00 
00 00 00 00 00 

where the origin is in the middle. We denote the maximum relative error for 
this class of neighborhoods by £^ where we restrict the limsup to vectors v with 
finite lengths W(v) (i.e. having coordiantes (x,y) with < y < x and p \ x). 
Our motivation for considering such neighborhoods is that it will turn out that 
(due to its special form) J\f c is easier to handle, but yields the same m.r. error 
as the corresponding neighborhood N c , in which 7V c (±p, 0) = N c (0, ±p) = p 
and N c (x,y) = cy/x 2 + y 2 otherwise ((x,y) G M p ). 

In Section 2 we introduce some notation and prove some preliminary results. 
In Sections 3 and 4 we compute the values of £^ and £® where £f is the 
maximum relative error £^ for optimal c and £^ = £ x p . We give all sequences 



lim sup 
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yielding minimal m.r. error in case of 5 x 5 and 7x7 neighborhoods, as well. 
In Section 4 we prove that Ep = £^ and = Ep and further show that 
Eg = Eg/ (2 + Eg) for all p. Finally, we draw some conclusions in Section 5. 



2 Definitions and basic properties 

Let N be a neighborhood defined on the mask M p . Put M* = M p \{(0, 0)}. We 
denote the value of N at position (n, k) by w(n, k) for (n, k) G M p . Throughout 
the paper we assume that w(±n,±k) = w(±k,±n) > for all (n,k) G M* 
and all possible sign choices. Hence it suffices to consider the values w(n,k) 
with < k < n < p. 

We can measure lengths of vectors and distances between points using neigh- 
borhood sequences. Note that such sequences provide a flexible and very use- 
ful tool in handling several problems in discrete geometry. For the basics and 
most important facts about such sequences, see e.g. the papers [20ll21f22f23ll24] 
and the references given there. Here we only give those notions which will be 
needed for our purposes. 

Let A = (iVj)?^ be a sequence of neighborhoods defined on M p and u, v G Z 2 . 
The sequence u = u , Ui, . . . , u m = v with — 6 M p is called an A- 

m 

path from u to v. The A-length of the path is defined as Wi(Hi — Ui-i). 

»=i 

The distance Wa(v — u) between u and v, which is the A- length of v — u, is 
defined as the minimal A-length taken over all A-paths from u to v. If the 
neighborhood sequence is fixed, then we suppress the letter A in the above 
notation. 

If Ni = N for all i, then the corresponding (constant) neighborhood sequence 
is denoted by A = N. We assume throughout the paper that for such sequences 
W(n,k) = w(n,k) holds for (n,k) G M p ; if it would not have been the case, 
then the function w := would have generated W, too. 

We call W a metric if for all u, v G Z 2 

• W{u) < oo {W is finite), 

• W{u) = <^ u = (W is positive definite), 

• W(u) = W(— u) {W is symmetric), 

• W(u + v) < W(u) + W{v) (W satisfies the triangle inequality). 

It follows from the above properties that W(u) > for every u G Z 2 . By our 
basic assumptions on w, every induced length function W is positive definite 
and symmetric. Furthermore, W satisfies the triangle inequality for u, v with 
u,v,u + v G M p by definition. 
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The first lemma shows that in case of a constant neighborhood sequence 
W(i/)/|?7| attains a minimal value which is reached already in M*. 

Lemma 1 Let N be a neighborhood defined on M p which induces the length 
function W on Z 2 . Then 

. . W(v) . w(v) 
hmmt — — — = mm . 

\v\-hx> \v\ vGM* \v\ 



PROOF. Let m = min ^ = ^ (« 6 M!). Then for all n we have 



V 



w ( n ^ = rn, so that liminf ^S^- < m. On the other hand, since > m for 

\nu\ ' \v\ — \v\ — 

every v G M*, it follows from the definition of shortest path and the triangle 
inequality for the Euclidean distance that 

W(v) > w(vi) = ^2 prf" ' Wi\ — 171 ^2 Kl — m \v\ 



for every v G Z 2 not equal to the origin. Thus lim inf > to. □ 

\v\^oo \ v \ 



The challenge is therefore to compute lim sup 



3 The maximum relative error for neighborhoods M c 



Let c be some positive real number with ^/ P 2+1 < c < 1. We shall study 

neighborhoods Af c on M p with Af c (n, k) — oo for which either n < p or k < 0, 
■A/^O 7 ) 0) = p and jV c (p, k) = c^/p 1 + /c 2 for < fc < p. We are interested in 
the length function W c induced by A c '■— 7\T C for points in the set {(x,y) G 
Z 2 : <y< x}. First we secure that under suitable conditions only two 
distinct steps occur in a shortest *4 c -path. 

Lemma 2 Let —j= — < c < 1. Then a shortest A c -path from (0, 0) to (mp, mr+ 
yp 2 +i 

k) with m,r, A; G Z, < r < p, < k < m consists only of steps (p, r) and 
(p,r + l). 



PROOF. Suppose a shortest path from (0, 0) to (mp, mr + k) with m,r,k G 
Z,,0 < r < p,0 < k < m contains two steps (p, t) and (p, -u) with t — 2 > -u > 0. 
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Replace the two steps with steps (p, t — 1) and (p, u + 1), and write Li and L 2 
for the length of the old and new paths, respectively. Then we have 



L l -L 2 > c^p 2 + t 2 - c^p 2 + (t- l) 2 + c^p 2 + u 2 - c^p 2 + (u + l) 2 



c(/ P (t) -f p (u + l)), 



where 



f p (x) = y/p 2 + x 2 




? + {x- l) 2 (z e Z> ). 



A simple calculation yields that is strictly monotone increasing in x, 

which shows that L\ — L 2 > 0. However, this contradicts the minimality of 
the length of the original path. 

Hence a shortest path may contain steps (p,t) and (p,t + 1) only, for some 
non-negative integer t. Since altogether we make m steps, this immediately 
gives that t = r, and our statement follows. □ 

Remark 3 The latter inequality is the most severe and explains why we re- 
strict c to values qreater than — . 

Corollary 4 Let Z — < c < 1 Then a shortest A c -path from (0, 0) to 
yp 2 +i 

(mp,mr) with < r < p consists of m steps (p,r). 

The next theorem gives the value of the approximation error for general p, in 
case of any neighborhood Af c on M p . 

Theorem 5 Let p > 1, — < c < 1. Then the m.r. error of A c to the 
Euclidean distance is given by 



PROOF. As a general remark we mention that to perform our calculations, 
we used the program package Maple ®0. 

Let p be a positive integer, and fix c with —t£ — < c < 1. As previously, it is 

Vp 2+1 

sufficient to consider the v4 c -length of points of the form (mp, k) where m is 
some positive integer and k is an integer with < k < mp. Write k = mq + r 
with < q < p and < r < m. The possible steps are (p, 0) of length p 
and (p,±i) of length Wi := C\fp 2 + i 2 (for |z| < p). From Lemma [2] and the 
inequalities p = W < W\ < . . . < W p we see that a path of minimal length 



Maple is a registered trademark of Waterloo Maple Inc. 



max(l — c, y 1 + c 2 + p 2 + c 2 p 2 — 2cpJp 2 + 1 — 1). 
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from (0,0) to a point (mp,mq + r) consists of r steps (p,q + 1) and m — r 
steps (p,q). Hence for the induced length function we get 

W(mp, mq + r) = rW q+ i + (m — r)W q . 



Put t = r/m, and recall that Wq = p and Wi = cy/p 2 + i 2 for % — 1, . . . ,p. Set 



_ ctVF+T + p(i-t) 



and for 1 < q < p 



Hg(t) = C-^ 1, 

/p 2 + (g + t) 2 



and let 

c ) = max |if,(*)| (0 < g < p) and /ip(p, c) = \H P (0)\. 

Now we investigate the error functions h q (p,c) for q = p, q = 0, < q < p, 
respectively. 

Suppose first that q = p. Then r = and k = mp. In this case we trivially 
have h p (p, c) = 1 — c. 

Assume next that q — 0. Then < k < p. Put 



t := p(cy / p 2 + 1 -p). 



A simple calculation yields that < t < 1, and that H is monotone in- 
creasing on the interval [0,io] an d monotone decreasing on the interval [to, 1]- 
Moreover, we have H o (0) = and H (l) = c — 1, hence H (t ) > 0. Thus we 
have 



h (p,c) = max(l-c, H (t )) = max(l-c, \Jl + c 2 +p 2 + c 2 p 2 - 2cp^p 2 + 1-1). 



Finally, suppose that < q < p, that is p < k < mp. Put 

t _ + yV + q 2 )(p 2 + (q +W) - p 2 - q 2 - q) 
{q + 1) ^r T ^_ q ^ p 2 + {q + 1) 2 

A simple calculation gives that < t q < 1, and that H q is monotone increasing 
on the interval [0, t q ], while monotone decreasing on the interval [t q , 1]. We also 
have H q (0) = H q {l) = c-l. Hence H q (t q ) < implies -H q (t q ) < 1 - c. Thus 
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we get 



h q (p,c) = max(l— c, H q (t q )) = max 



1 — c, c 



V 



\ 



- 1 



l + Wl- 



Now we calculate the error function 



h(p, c) := limsup 

p\n, n>k>0 



Vn 2 + A; 2 



max c). 



0<9<P 



Observe first that for fixed p and c the function c) is monotone decreasing 

in q with 2 < q < p. Hence c) < /ii(p, c) for q — 2, . . . ,p. Further, again 

by Maple, we obtain that for any c with , p < c < 1 

Vp + 1 



1 



1 + a/1 



< a/1 + c 2 + p 2 + c 2 p 2 - 2cpJp 2 + 1 



(p 2 + l)(p 2 +4) 

holds, which implies hi(p,c) < h (p,c). Hence 



h(p, c) = max(l — c, \j\ + c 2 + p 2 + c 2 p 2 — 2cp\Jp 2 + 1 — 1) 



and the theorem follows. □ 



The following corollaries provide the m.r. errors (when c 
(when c = 1), respectively. 

Corollary 6 Let p be a positive integer. Then we have 



Cp) and £p 



D 



pVp^Tl + 2^p 2 + 1 - pvV + 1 - 2 

^2 



TTierf is, £/ie sequence A = A c b of period p given by A = N c b yields the smallest 
m.r. error among all sequences A c of period p. Moreover, the error is given by 



s b = 1 _ c b = p 2 + 2 - pVF+t - 2 yp 2 + 1 - pVp^+t 



V 



0.0858 



L5 ~^ i cfM.--- 

J9 2 VP 4 / 
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PROOF. Put 



f(c) = 1 — c and g(c) = J 1 + c 2 + p 2 + c 2 p 2 — 2cp\jp 2 + 1 — 1. 

A straightforward computation shows that / is strictly monotone decreasing, 

while g is strictly monotone increasing for X — < c < 1. Hence there is a 

vV+i 

unique solution of the equation /(c) = g(c) in this interval. By Theorem [5] 
this solution is given by 

B _ pVp^Tl + 2^ 2 ^V^p~^ 2= T\ - 2 
Thus the statement follows. □ 

Corollary 7 Let p be a positive integer. Then the sequence A = A± of period 
p given by A = ~Ni (corresponding to the choice c = 1 ) has m.r. error 

PROOF. On substituting c = 1 into the formula of Theorem[5], the statement 
follows immediately. □ 



Now we give the best approximating sequences realizing the minimal maxi- 
mum relative error for 5x5 matrices (p = 2) in Theorem [8] and for 7x7 
matrices (p = 3) in Theorem IHJ respectively. 

Theorem 8 Let < c < 1. Let A c = ~7J~ C be the corresponding sequence on 
Mi- Then the minimal m.r. error to the Euclidean distance among the neigh- 
borhood sequences A c is attained if and only if 



-2 ) 



W\ = s and u < W 2 < v, 



where 



5-^+^25-10^ 

s = i — 2.1943, 

2 

2^2 2s 
u = — ^s R3 2.7756 and v = 2 + — ^ 2.8777. 
y/E 5 



Further, the m.r. error is given by 



D n S 3-^-^5-2^ 

p b = 1 - c f = 1 - = V — w 0.0187. 

2 2 2 
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PROOF. For any even n with < k < n the possible steps are (2, 0) of 
length 2, (2, 1) and (2, -1) of length W x , and (2, 2) and (2, -2) of length W 2 . 
From Lemma [2] and the inequality 2 < W\ < W 2 we see that the path from 
(0, 0) to (n, k) of minimal length consists of k steps (2, 1) and § — k steps (2, 0) 
if < k < n/2 and of k - n/2 steps (2, 2) and n-k steps (2, 1) if § < k < n. 
Hence we have for the induced length function 



W(n,k) 



kWx + n- 2k, if k < § , 

(n - k)Wi + (k — | )W 2 , otherwise. 

Put t = k/n. Then the error function is given by 

W(n, k) 



h(Wi, W 2 ) := limsup 

2|n, n>fe>0 



max max 

\0<i<± 



t{W 1 - 2) + 1 



max 

i ^ +^ 1 



y/n 2 + k 2 
1 - 1) W 1 + {t - \)W 2 



l 



Our aim is to choose W\ and W 2 such that h(W\, W 2 ) is minimal. For fixed 
Wi, define the function H : R> R by 



#o(*) 



- 2) + 1 



Put t = Wi — 2. We observe that H is monotone increasing on [0,to] and 
monotone decreasing on [to, oo). Hence, as Hq(0) = 1, 



max(\H (t) - 1|) = max (ifofo) - 1, 1 - H ( - 



o<t<| 



max - 4Wi + 5 - 1, 1 - ^| 



if < 5/2 and 



max(\H (t) - 1|) = # Q - 1 = -± - 1 



otherwise. Clearly, 



min (/i(Wi, W 2 )) > minmax WW-, 2 - 4PFi + 5-1, 



1 - 



Wi 

7H 



(i) 



A calculation gives that the minimum of the right-hand side is achieved for 



5-^+1/25-10^ 
Wi = s:= ^ — « 2.1943 
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and equals 



Vs 2 - 4s + 5 - 1 = 1 - — 

VS 



VS-J5-2VS 



0.0187. 



Now we fix the value s of W±, and show that we can choose W 2 in a way to have 
equality in ([[]). In fact we completely describe the set of the appropriate W2-S. 
Consider the maximum over t G [1/2,1]. For fixed W 2 , define the function 
H x : R> ->■ R by 

(l-t)Wi + (*-§)W 2 



#i(*) 



Observe that H\ attains its maximum at ti := 2 2 ^ 1 ZyyJ (which is positive) 
and further, H\ is monotone increasing in [0, ti] and monotone decreasing in 
[ti, 00). Hence 



max (\Hi(t) - 1|) = max 1 - H x 

i<t<i 



Hi{ti 



l,l-^i(l 

Wi J(2Wi-w 2 y + A(w 2 -Wiy w 2 

max 1 =. 1,1 

VS 2 ' 2^2, 

if 1/2 < ti < 1, and 

max (\Hi(t) - 1|) = max (\H X (\\ - 1|, |^(1) - 1 



max 1 — 



Wx 

VS' 



2V2 

otherwise. By our choice of W\, we have that 



Wo 



1 - 



7E 



1 - 



vs 



0.0187. 



The values of 



w 2 _ 1 
2V2 1 



and 



y /(2W 1 -W 2 ) 2 +4,(W 2 -Wtf 



do not exceed this 



value if and only if u < W 2 < v where u and v are defined in the statement 
of the theorem. We conclude that h(Wi,W 2 ) attains its minimum 1 — if 
Wi — s and u <W 2 <v. 

The above argument shows that £f = 1 — Hence the minimum among 
neighborhoods J\f c is realized for c = cf = ^= and for no other value of c. □ 

Theorem 9 Let < c < 1. Le£ *4. c = A/^ fre the corresponding sequence on 
M3. Then the minimal m.r. error to the Euclidean distance among the neigh- 
borhood sequences A c is attained if and only if 



c = cf , Wx = s, u < W 2 < v, q < W 3 < r, 



12 



where 



30 - 2^\Q + 2^100 - 30^ 

s = £ — « 3.1340, 

9 



13 

u = -, — s ~ 3.5733, 
10 



143s - 18V13 + 6V1690 - 143V13s 

v = ¥ « 3.5944, 

121 

3s 

g = -j= « 4.2047, 
\/5 



3v/l3s 2 -52v / 10s + 520 - 10PV| 15 W 2 
r = — p= h 



13^10 13 ' 

and in the definition of r, W 2 can be any number with u < W 2 < v. Further, 
the m.r. error is given by 



r> r> s 11-3-n/10-2V10-3-n/10 
£f = l- c? = l- -== = y — « 0.0089. 

3 3 9 



PROOF. Let 3\n and < k < n. The possible steps are (3,0) of length 3, 
(3, ±1) of length W u (3, ±2) of length W 2 , and (3, ±3) of length W 3 . From the 
inequalities 3 < W\ < W 2 < W 3 it follows that the path from (0, 0) to (n, k) 
of minimal length consists of k steps (3, 1) and | — k steps (3, 0) if < k < |; 
of k - § steps (3, 2) and f - fc steps (3, 1) if § < k < f ; of k - f steps (3, 3) 
and n — k steps (3, 2) if < k < n. Hence we have for the induced length 
function 



W(n, k) 



kW 1 + n-3k, if k < n/3, 

(2n/3 - k)Wi + (k - n/3)W 2 , if n/3 < k < 2n/3, 

(n - k)W 2 + (k- 2n/3)W 3 , otherwise. 



Put t = k/n, and define the functions Hi : IR> — > K (i = 0, 1, 2) by 
and 

(l-t)H/ 2 + (t-|)H/ 3 

^ = 7TT¥ • 

Then for fixed W 1: W 2: W 3 the error of approximation is given by 
h{W u W 2 ,W 3 ) — max I max \Ho(t) — 1|, max \H\[t) — 1|, max \H 2 (t) — 1| ] . 

Vo<t<i i<t<§ l<t<l J 
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Let 



t = W 1 - 3, h 



3(W 2 -W 1 ) + 3(W 3 -W 2 ) 



2W X -W 2 3W 2 -2W 3 
and observe that all to, ti and t 2 are positive. By differentiation and following 
standard calculus, we get that for i = 0,1, 2, Hi is monotone decreasing if ti £" 
[i/3, (i + l)/3], and that Hi is monotone increasing in [i/3, U] and monotone 
decreasing in [ti, (i + l)/3] otherwise. Hence from H (0) = 1 we get that 



max (\H (t) — 1|) = max H (to 



0<t< 



II ( - 



1 



max 



Hence obviously, 



W^-6W 1 + 10-1, 



- 1 



min h{W x , W 2 , W 3 ) > minmax JW? - QW Y + 10-1, 

Wl,W2,W3 Wl \ v 



1 - 



Wi 



(2) 



By a simple calculation we get that the minimum of the right-hand side is 
achieved for 



30-2V10 + 2V100 - 30V10 

W x = s:= y — w 3.1340 

9 



and equals 



M 



Vs 2 - 6s + 10 - 1 = 1 



3^+7 + 2^10-3^ 
9^ 



0.0089. 



Now we fix the value s of W\, and show that we can choose W 2 and W3 in a way 
to have equality in (j2J). More precisely, we completely describe the set of the 
appropriate pairs (W 2 , W3). For this purpose, first we consider the maximum 
of H\ over t £ [1/3,2/3]. In a similar manner as in the proof of Theorem [HI 
we obtain that 



max (\Hi(t) — 1|) = max 



3 L <*<"3 



1 



max 



- 1 



(2Wi - Wa) 2 + 9(W 2 - Wi, 



-1, 



- 1 



Using our choice for Wi, a simple calculation gives that the above maximum 
does not exceed the value of M precisely when u < W 2 < v, where u and v 
are defined in the statement of the theorem. So let W 2 be any fixed number 
from the interval [u,v], and consider the the maximum of H 2 over t £ [2/3, 1]. 
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Now we get that 

max (\H 2 (t) — 1|) = max I 

'i<t<i V 



Ho 



,H 2 (t 2 )-l,\H 2 (l)-l\) = 



max 



Wo 



13 



- 1 



(3W 2 - 2W 3 ) 2 + 9(W 3 - W 2 f 



1, 



3V% 



- 1 



Using our choice for W 1 and W 2 , a simple calculation yields that the above 
maximum is not larger than M if and only if q < W 3 < r, where q and r are 
given in the statement. (Note that 4.2766 < r < 4.2804.) 

The above argument shows that £f = 1 — Hence the minimum among 
neighborhoods Af c is realized for c = cf = and the theorem follows. □ 



4 Equivalence of m.r. errors for M p neighborhoods 



In this section we compute the m.r. errors Ep, Ep and E°. First we introduce 
neighborhoods N c on M p defined by iV c (0,0) = oo, N c (n, 0) = iV c (0,n) = \n\ 
for < \n\ < p, N c (n,k) = c^n 2 + k 2 for (n, k) G M p ,nk ^ 0. Let W c 
denote the length function induced by the sequence N c . We show that the 
corresponding m.r. error E c satisfies E c = S c for every considered value of c. It 
then follows that Ep 3 = E p and E® = E p for every p > 1. 

Lemma 10 Let X — < c < 1. There is a shortest N c -path from (0,0) to 
Vp 2 + 1 

(mp, k) with < k < m which consists of steps of the form (p, 0) and (p, 1). 



PROOF. Suppose a shortest path from (0, 0) to (mp, k) contains a step (g, h) 
with h < 0. Then it also contains a step (i, j) with j > 1. But it is shorter to 
replace both steps with steps (g, h+1) and (i, j — 1). A similar argument can 
be used to exclude steps (g, h) with h > 1. So every shortest path from (0, 0) 
to (mp, k) contains only steps of the forms (g, 0) and (g, 1). 

If k = m, then taking only steps (p, 1) gives the shortest path length because 
of the triangle inequality for the Euclidean distance and the inequality c < 1. 
Suppose that there is a step (g, 1) with g < p in a shortest path from (0, 0) to 
(mp, k) with < k < m. Then there is also a step (h, 0) with h > 0. But we 
can replace both steps with steps (g + 1, 1) and (h — 1, 0) and make the path 
shorter. Therefore all the steps of the form (g, 1) are of the form (p, 1). The 
remaining steps can be combined to steps of the form (p, 0). □ 
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Lemma 11 Let p be fixed. Let — < c < 1. The m.r. error of the neigh- 

a/p 2 +i 

borhood sequence N c is equal to if c = 1 and equal to S,^ if c assumes the 
value Cp from Corollary® 



PROOF. Because of symmetry it suffices only to consider points (n, k) with 
< k < n. First let c = 1. By definition N(n, k) = \/n 2 + k 2 for (n, k) G M p . 
Hence the induced length function satisfies Wi(n, k) > \(n, k)\ for all (n, k) e 
Z 2 . Thus 

mm = > 1 

\/n 2 + k 2 

where the minimum is taken over all (n, k) e Z 2 with (n, k) ^ (0, 0). On the 
other hand, by Lemma HO, the shortest Ni path from (0, 0) to (mp, k) with 

< k < m consists of steps of the forms (p, 0) and (p, 1) which have lengths 
p and i/p 2 + 1, respectively. Hence Wi(mp,k) = Wi(mp,k) for < k < m. 
If n = mp + r with < r < p, then |Wi(n, k) — Wi(mp, k)\ < p. Note that 
in view the proof of Theorem [5] (in particular, since h^(p,c) > hi(p,c) for all 

1 > 1 there) we have 

Wi(mp, k) W 1 (mp,k) 
hmsup — = hmsup — — . 

|(m P ,fc)Koo \{mp,k)\ \{ mp ,k)\-+<x> \{mp,k)\ 

0<k<mp 0<k<m 

Thus on the one hand it follows that 

Wifak) Wiimp^k) 
hmsup — — — = hmsup — — > 

|(n,fc)|->oo \{n,k)\ |(mp,fc)[-+t» \{mp,k)\ 

0<k<n 0<k<mp 

^ Wi(mp,k) ,. Wifmp, fc) .. Wi(mp, fc) 
> hmsup — — = hmsup — — = hmsup — — . 

\(mp,k)\->oc \{mp,k)\ \(mp,k)\-xx> \{ m P> k )\ \(mp,k)\-xx> \( m Pi k )\ 
0<k<m 0<k<m 0<k<mp 

On the other hand, by W\(mp, k) < Wi(mp, k) for all m,p and k, we also 
have that 

Wi(mp,k) , WUrnp, k) WAn^k) 
hmsup — — — > hmsup — — — = hmsup 



|(mp,fc)Koo \(mp,k)\ |( m p,fc)|->oo \(mp,k)\ |(n,fc)|->oo \(n,k)\ 

0<k<mp 0<k<mp 0<k<n 



Hence 



and by 



Wi(n,k) .. Wi(mp,k) 
hmsup — — — = hmsup — — 

\(n t k)\-xx> \( n ' k )\ \{mp,k)\-¥oo \{mp,k)\ 
0<k<n 0<k<mp 



limSUp ^TT—^JTT =l+£p, 

\(mp,k)\-+ao \[mp,k)\ 

the m.r. error of Ni equals 5r . 
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Next let c = = 1 — £f . Then — < c < 1, and, by construction, 

Wc{p, 0) = p, W c (p, k) = cy/p' 2 + k 2 for < k < p, and W c (n, k) = c\fn 2 + k 2 
for < k < n < p. Hence 

• W c (n, k) B 
mm . = = c = 1 — t„ . 

(n,k)eM; ^n 2 + k 2 p 

Thus 

W c b(u, k) 
li m j n f p = i _ £ B 

|(n,fc)Koo \(n,k)\ p ' 

On the other hand, by Lemma [TDl the shortest iV c path from (0, 0) to (mp, k) 
with < k < m consists of steps of the form (p,0) and (p, 1). By a similar 
reasoning as above we obtain that 

W c s(n,k) 
hmsup " =1 + 8 

|(n,fc)|->oo \{ n ,k)\ 

Thus the m.r. error of N c equals Sp . □ 

Theorem 12 For every p we have EL 3 = Si 3 and EL = . 



PROOF. We first consider the D-c&se. Suppose the neighborhood iV on M p 
induces a length function W : Z 2 — > R>o such that W(v) > \v\ for all i?£Z 2 
and W has m.r. error i?^ 3 . It can only improve the m.r. error if we replace the 
value N(n, k) for some (n, k) e M* with a smaller value > |(n, k)\. Therefore 
we may assume without loss of generality that N = N±. Hence E® = Sp. 

Now we turn to the 5-case. Suppose a neighborhood iV on M p induces a length 
function W such that W(n, 0) = W(0,n) = \n | for n e Z and (1 - E^)\v\ < 
W(v) < (1 + .Ejf for all -u G 1? . Without loss of generality we may replace 
all values N(n,k) for (n,k) G M* with \n\ if = 0, with \k\ if n = 0, 
and with (1 — Ep)\(n,k)\ otherwise. Thus Ei 3 equals the m.r. error of the 
neighborhood sequence Ni_ e b. We know from Lemma [11] and Corollary [6] 
that if c = Cp , then the m.r. error of iV c equals 8^ = 1 — cp 3 . Hence Ep 3 < Sp 3 . 
From N(n, k) > (1 - Ef)\(n,h)\ > cg\(n,k)\ for all (n,k) G M* we obtain 
W(v) > W c b(v) for all v G Z 2 . Hence 

R W(v) W c b(v) 



1 + E„ = inf lim sup — --— > lim sup — ^- — = 1 + 8, 

" \v\-+oo \ u \ \v\-*ao 



by Lemma EJ Thus Ep =8^. □ 

Finally, we compute the minimal m.r. error Ep for the class of arbitrary neigh- 
borhoods N defined on M p . Observe that the m.r. error Ep is attained by the 
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length function W corresponding to the neighborhood N defined by w(v) = 
(1 — Ep)\v\ for v £ M*, since should not assume a smaller value than 
1 — Ep and the limsup- value cannot increase if we decrease some w(v). Clearly, 
the length function W corresponding to N is just \ c W\ where W\ is the 

length function on N\. Recall that N\ has m.r. error . Therefore we have 
1 + E° = hmsup ^ = (1 + <)(1 - E°). (3) 

\v\— >oo \ v \ 

By a simple calculation we get E^ = 2+ v e d ■ So we have proved 
Theorem 13 For every p > 1 we have 

E c E p _ y/2p 2 + 2 - 2 Pv ^ + T - 1 l |g 
P 2 + E p J2p 2 + 2-2p^p r TT+l 1Q P 2 




Remark 14 Observe that E^ is about 37% larger than E^. This is the price 
to be paid for the restriction W(n, 0) = \n\ for n e Z. T/ie fa/ne of E^ is 
about twice the error E^ . This is due to the fact that the negative and positive 
deviations in E^ are added to the positive deviation in E® . 



5 Conclusion 



In this paper, we have determined the smallest possible maximum relative 
error of chamfer distances with respect to the Euclidean distance under vari- 
ous conditions. We have dealt with approximating distances from three main 
aspects: supposing that a horizontal/vertical step has a weight 1 in the local 
chamfer neighborhoods, majorating the Euclidean distance, and also without 
any constraint. We have calculated optimal weights for small (5x5 and 7x7) 
neighborhoods in a certain case, as well. Our framework is embedded in the 
theory of neighborhood sequences with possible generalizations in this field. 
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